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Abstract 

The Baikal Neutrino Telescope NT200 takes data since April 1998. On April 9th, 2005, the 10 Mton scale detector 
NT200+ was put into operation in Lake Baikal. Selected results obtained during 1998-2002 with the neutrino 
telescope NT200 are presented. 
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1. Introduction 

The Baikal Neutrino Telescope NT200 takes 
data since April 1998. On April 9th, 2005, the 10- 
Mton scale detector NT200+ was put into opera- 
tion in Lake Baikal. Description of site properties, 
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detector configuration and performance have been 
described elsewhere [1,2,3,4]. 

In this paper we present main results obtained 
from the analysis of data taken with the Baikal 
neutrino telescope NT200 between April 1998 and 
February 2003. 
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2. Selected results obtained with NT200 

2.1. Atmospheric neutrinos 

The signature of charged current muon neutrino 
events is a muon crossing the detector from below. 
Muon track reconstruction algorithms and back- 
ground rejection have been described elsewhere [5] . 
Compared to [5] the analysis of the 4-year sample 
(1038 days live time) was optimized for higher sig- 
nal passing rate, and accepting a slightly higher 
contamination of 15-20% fake events [6]. A total 
of 372 upward going neutrino candidates were se- 
lected. From Monte-Carlo simulation a total of 385 
atmospheric neutrino and background events are 
expected. The skyplot of these events is shown in 
Fig. 1. 
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Fig. 1. Skyplot (galactic coordinates) of neutrino events 
for five years. The solid curve shows the equator. 

2.2. Search for Neutrinos from WIMP 
Annihilation 

The search for WIMPs with the Baikal neutrino 
telescope is based on a possible signal of nearly 
vertically upward going muons, exceeding the flux 
of atmospheric neutrinos. The method of event se- 
lection relies on the application of a series of cuts 
which are tailored to the response of the telescope 
to nearly vertically upward moving muons [7] . The 
applied cuts select muons with -1< cos(#) <-0.75 
and result in a detection area of about 1800 m 2 
for vertically upward going muons. The energy 
threshold for this analysis is -Ethr ~ 10 GeV i.e. 
significantly lower then for the analysis described 
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Fig. 2. Angular distributions of selected neutrino candi- 
dates as well as expected distributions in a case with and 
without oscillations (solid and dashed curves respectively) . 

in section 2.1 (E thr ~ 15-20 GeV). Therefore 
the effect of oscillations is stronger visible. We 
expect a muon event suppression of (25-30)% due 
to neutrino oscillations assuming dm 2 =2.5T0 -3 
eV 2 with full mixing, 6 rn w 7r/4. 

From 502 days of effective data taking between 
April 1998 and February 2000, 24 events with -1< 
cos(#) <-0.75 have been selected as clear neutrino 
events. The angular distribution of these events 
as well as the MC - predicted distributions are 
shown in Fig. 2. For the MC simulations we used 
the Bartol96 atmospheric neutrino flux [8] without 
(dashed curve) and with (solid curve) oscillations. 
Within lcr statistical uncertainties the experimen- 
tal angular distribution is consistent with the pre- 
diction including neutrino oscillations. 

Regarding the 24 detected events as being in- 
duced by atmospheric neutrinos, one can derive an 
upper limit on the additional flux of muons from 
the center of the Earth due to annihilation of neu- 
tralinos - the favored candidate for cold dark mat- 
ter. The 90% C.L. muon flux limits for six cones 
around the opposite zenith as well as muon flux 
limits for different neutralino masses obtained with 
NT200 (E thl >10 GeV) in 1998/99 are shown in 
Fig. 3 and Fig. 4, and compared to limits obtained 
by Baksan [9], MACRO [10], Super-Kamiokande 
[11] and AMANDA (from the hard neutralino an- 
nihilation channels) [12]. 
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3. Limits on the excess muon flux from the center of 
Earth versus half-cone of the search angle. 



Nhit > 30. In order to reduce the background from 
downward atmospheric muons we restrict ourself 
to monopoles coming from the lower hemisphere. 
For an upward going particle the times of hit chan- 
nels increase with rising z-coordinates from bot- 
tom to top of the detector. To suppress downward 
moving particles, a cut on the value of the time— z- 
correlation, Ct z , is applied: 



° tz ~~ AT > U 



(1) 



where i, and Z{ are time and z-coordinate of a fired 
channel, t and z are mean values for times and z- 
coordinates of the event and at and a z the rms- 
errors for time and z-coordinates. 
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Fig. 4. Limits on the excess muon flux from the center of 
the Earth as a function of WIMP mass. 

2.3. A search for fast magnetic monopoles 

Fast magnetic monopoles with Dirac charge g = 
68. 5e are interesting objects to search for with deep 
underwater neutrino telescopes. The intensity of 
monopole Cherenkov radiation is w 8300 times 
higher than that of muons. Optical modules of the 
Baikal experiment can detect such an object from 
a distance up to hundred meters. The processing 
chain for fast monopoles starts with the selection 
of events with a high multiplicity of hit channels: 
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Fig. 5. Ctz distributions for experimental events (tri- 
angles), simulated atmospheric muon events (solid), and 
simulated upward moving relativistic magnetic monopoles 
(dotted); multiplicity cut N^u > 30. 



In Fig. 5 we compare the Ct z -distribution for 
experimental data (triangles), simulated atmo- 
spheric muon events (solid curve) with simulated 
upward moving monopole events (dotted curve). 

Within 994 days of live time using in this anal- 
ysis, about 3 ■ 10 s events with Nhit > 4 have been 
recorded, with 20943 of them satisfying cut 
(Nhit > 30 and Ctz > 0). For further background 
suppression (see [13] for details of the analysis) we 
use additional cuts, which essentially reject muon 
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events and at the same time only slightly reduce 
the effective area for relativistic monopoles 1 : 

(i) N hit > 35 and C tz > 0.4 + 0.6 

(ii) x 2 determined from reconstruction has to be 
smaller than 3 

(iii) Reconstructed zenith angle 9 > 100 deg 

(iv) Reconstructed track distance from NT200 
center R rec > 20 -f- 25 m. 

No events from the experimental sample pass 
cuts 1-4. The acceptances A e f f for monopoles with 
(3 =1, 0.9 and 0.8 have been calculated for all 
NT200 operation configurations (various sets of 
operating channels) . For the time periods included, 
A e f f varies between 3 TO 8 and 6 • 10 8 cm 2 sr (for (3 — 
1). From the non-observation of candidate events 
in NT200 and the earlier stage telescopes NT36 
and NT96 [14], a combined upper limit on the flux 
of fast monopoles with 90% C.L. is obtained. 

In Fig. 6 we compare this upper limit for an 
isotropic flux of fast monopoles obtained with the 
Baikal neutrino telescope to the limits from the un- 
derground experiments Ohya [15] and MACRO[16] 
and to the limit reported for the underice detector 
AMANDA [17]. The Baikal limit is currently the 
most stringent one. 
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Fig. 6. Upper limits on the flux of fast monopoles obtained 
in this analysis (Baikal) and in other experiments. 



1 Different values of cuts correspond to different NT200 
operation configurations. 



2.4. A search for extraterrestrial high-energy 
neutrinos 

The BAIKAL survey for high energy neutrinos 
searches for bright cascades produced at the neu- 
trino interaction vertex in a large volume around 
the neutrino telescope [2]. We select events with 
high multiplicity of hit channels N^t, correspond- 
ing to bright cascades. To separate high-energy 
neutrino events from background events a cut to 
select events with upward moving light signals has 
been developed. We define for each event t m i n = 
min(£i — tj), where ti, tj are the arrival times at 
channels i, j on each string, and the minimum over 
all strings is calculated. Positive and negative val- 
ues of t m [ n correspond to upward and downward 
propagation of light, respectively. 

Within the 1038 days of the detector live time 
between April 1998 and February 2003, 3.45 x 
10 s events with N^ t > 4 have been recorded. 
For this analysis we used 22597 events with hit 
channel multiplicity N^ it >15 and i mm >-10 ns. 
We conclude that data are consistent with simu- 
lated background for both £ m ; n and N^ it distribu- 
tions. No statistically significant excess above the 
background from atmospheric muons has been ob- 
served. To maximize the sensitivity to a neutrino 
signal we introduce a cut in the (t m [ n , N^ it ) phase 
space. 

Table 1 

Expected number of events A r mo( j e [ and model rejection 
factors for models of astrophysical neutrino sources. 



BAIKAL 


AMANDA [18,19,20] 


Model u e 


+ v„ 4 


■v T "90%/^modcl 




1CT 6 x E~ 2 


3.08 


0.81 


0.86 


SS Quasar [21] 


10.00 


0.25 


0.21 


SP u [22] 


40.18 


0.062 


0.054 


SP 1 [22] 


6.75 


0.37 


0.28 


P PI [23] 


2.19 


1.14 


1.99 


M pp + P7 [24] 


0.86 


2.86 


1.19 


MPR [25] 


0.63 


4.0 


4.41 


SeSi [26] 


1.18 


2.12 





Since no events have been observed which pass 
the final cuts upper limits on the diffuse flux of 
extraterrestrial neutrinos are calculated. For a 90% 
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Fig. 7. Left panel: all-flavor neutrino flux predictions in different models of neutrino sources compared to experimental upper 
limits to E~ 2 fluxes obtained by this analysis and other experiments (see text). Also shown is the sensitivity expected for 3 
live years of the new telescope NT200+ [4,34]. Right panel: Baikal experimental limits compared to two model predictions. 
Dotted curves: predictions from model SS [21], SeSi [26] and SS05 [35]. Full curves: upper limits to spectra of the same 
shape. Model SS is excluded (MRF=0.25), model SeSi is not (MRF=2.12). 



confidence level an upper limit on the number of 
signal events of ng % =2.5 is obtained assuming 
an uncertainty in signal detection of 24% and a 
background of zero events. 

A model of astrophysical neutrino sources, for 
which the total number of expected events, iV mo( j e i , 
is large than n 90 %, is ruled out at 90% CL. Table 1 
represents event rates and model rejection factors 
(MRP) n 90 %/N mo d e \ for models of astrophysical 
neutrino sources obtained from our search, as well 
as model rejection factors obtained recently by the 
AMANDA collaboration [18,19,20]. 

For an E~ 2 behaviour of the neutrino spectrum 
and a flavor ratio v & : : v T = 1 : 1 : 1, the 90% 
C.L. upper limit on the neutrino flux of all flavors 
obtained with the Baikal neutrino telescope NT200 
(1038 days) is: 

E 2 $ < 8.1 x 10- 7 cm" 2 s" 1 sr" 1 GeV. (2) 

For the resonant process with the resonant 
neutrino energy Eq — 6.3 x 10 6 GeV the model- 
independent limit on v e is: 



< 3.3 x lO^curW^GeV- 1 . (3) 

Fig. 7 (left panel) shows our upper limit on the 
all flavor E~ 2 diffuse flux (2) as well as the model 
independent limit on the resonant v e flux (dia- 
mond) (3). Also shown are the limits obtained by 
AMANDA [18,19,20] and MACRO [27] , theoretical 
bounds obtained by Berezinsky (model indepen- 
dent (B) [28] and for an E~ 2 shape of the neutrino 
spectrum (B(E~ 2 )) [29], by Waxman and Bahcall 
(WB) [30], by Mannheim et al.(MPR) [25], predic- 
tions for neutrino fluxes from topological defects 
(TD) [26], prediction on diffuse flux from AGNs 
according to Nellen et al. (NMB) [31], as well as 
the atmospheric conventional neutrino fluxes [32] 
from horizontal and vertical directions ( {y) upper 
and lower curves, respectively) and atmospheric 
prompt neutrino fluxes {v pr ) obtained by Volkova 
et al. [33] . The right panel of Fig. 7 shows our upper 
limits (solid curves) on diffuse fluxes from AGNs 
shaped according to the model of Stecker and Sala- 
mon (SS, SS05) [21,35] and of Semikoz and Sigl 
(SeSi) [26], according to Table 1. 
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3. Conclusion 

The Baikal neutrino telescope NT200 is taking 
data since April 1998. The upper limit obtained 
for a diffuse (v e + + v T ) flux with E~ 2 shape is 
E 2 <& = 8.1 x 10 -7 cnr~ 2 s _1 sr _1 GeV. The limits on 
fast magnetic monopoles and on a v e flux at the res- 
onant energy 6.3xl0 6 GeV are presently the most 
stringent. To extend the search for diffuse extrater- 
restrial neutrinos with higher sensitivity, NT200 
was significantly upgraded to NT200+, a detec- 
tor with about 5 Mton enclosed volume, which 
takes data since April 2005 [4,34]. The three-year 
sensitivity of NT200+ to the all-flavor neutrino 
flux is approximately 2 x 10 _7 cm~ 2 s _1 sr _1 GeV 
for E >10 2 TeV (shown in Fig. 7). In parallel with 
exploiting NT200+ we started research & develop- 
ment activities towards a Gigaton Volume Detec- 
tor in Lake Baikal. 
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